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We compute the pion light-cone wave function and the pion quark distribution amplitude in 
the Nambu- Jona-Lasinio model. We use the Pauli-Villars regularization method and as a result 
the distribution amplitude satisfies proper normalization and crossing properties. In the chiral 
limit we obtain the simple results, namely tp-nix) = 1 for the pion distribution amplitude, and 
J d 2 k± } & 7r (x,k±)k 2 L = (k±) = —M(uu)/f% for the second moment of the pion light-cone wave 
function, where M is the constituent quark mass and / T is the pion decay constant. After the QCD 
Gegenbauer evolution of the pion distribution amplitude good end-point behavior is recovered, and 
a satisfactory agreement with the analysis of the experimental data from CLEO is achieved. This 
allows us to determine the momentum scale corresponding to our model calculation, which is close to 
the value Qo = 313 MeV obtained earlier from the analogous analysis of the pion parton distribution 
function. The value of (k±) is, after the QCD evolution, around (400 MeV) 2 . In addition, the 
model predicts a linear integral relation between the pion distribution amplitude and the parton 
distribution function of the pion, which holds at the leading-order QCD evolution. 
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I. INTRODUCTION 



The study of high-energy exclusive processes [|l] pro- 
vides a convenient tool of learning about the quark sub- 
structure of hadrons. In this limit the total amplitude 
factorizes into a hard contribution, computable from per- 
turbative QCD, and a soft matrix element which requires 
a non-perturbative treatment. From the point of view of 
chiral symmetry breaking a particularly interesting pro- 
cess is provided by the 7* — > 7*7r° transition form fac- 
tor. For real photons its normalization is fixed by the 
anomalous breaking of chiral symmetry by the tt° — > 77 
decay. In addition, in the limit of large photon virtual- 
ities, factorization allows us to define the leading-twist 
pion distribution amplitude as a low energy matrix el- 
ement whose normalization is fixed by the pion weak- 
decay constant, a spontaneous chiral symmetry break- 
ing feature of the QCD vacuum. It seems obvious that 
such a process offers a unique opportunity not only to 
learn about the interplay between high and low ener- 
gies, but also to understand the relation between the 
spontaneous and the anomalous chiral symmetry break- 
ing. Radiative logarithmic corrections to the pion dis- 
tribution amplitude (PDA) can be easily implemented 
through the QCD evolution equations [||, ||, which yield 
for Q 2 — > 00 the asymptotic wave function of the form 
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tpn(x, 00) = 6x(l — x). Moreover, the pion transition 
form factor has been measur ed by the CELLO @ and, 
recently, the CLEO collaborations Q . A theoretical anal- 
ysis of PDA based on these data and light-cone sum rules 
has been undertaken ||, showing that at Q — 2.4 GeV 
PDA is neither asymptotic, nor possesses the double- 
hump structure |7| proposed in early works || fi~9fl . 

The pion distribution amplitude has been evaluated 
with QCD sum rules f§, [lOj, |l[ |§, || @| 5 in stan- 
dard 1 15 1 (only the second f-moment) and transverse lat- 
tice approaches Jit] , [l~7|. [l8| , and in chiral quark mod- 

els g, || g^gfgTll §,§§!■ ^ chi- 
ral quark models the results are not always compatible 
to each other, and even their interpretation has not al- 
ways been the same. While in same cases there are 
problems with chiral symmetry and pro per normaliza- 
tion ||, |5 §5), m other cases [|§ ||, g, || |§ |§ || 
it is not clear how to associate the scale at which the 
model is defined, necessary to define the starting point for 
the QCD evolution. Nevertheless, there is a precise way 
to identify the low energy scale, Qo, at which the model is 
defined, namely the one at which the quarks carry 100% 
of the total momentum |30| . |3l| ■ The fact that several 
calculations §0|, |l[ [|§ [24[ ||7|^, || produce a PDA 
strongly resembling the asymptotic form suggests that 
their working scale is already large, and the subsequent 
QCD evolution becomes unnecessary, or numerically in- 
significant. This also tacitly assumes that these models 
already incorporate the QCD radiative corrections. 

In the present paper we compute the pion distribution 
amplitude and the pion light-cone wave function within 
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the Nambu-Jona-Lasinio (NJL) model |32, |33| in a semi- 
bosonized form using the Pauli-Villars (PV) regulariza- 
tion method j34j|. This method has been introduced in 
Refs. |p5 , 36 in the context of chiral perturbation theory, 
as well as for chiral solitons. From the point of view of 
the NJL model the study of exclusive processes becomes 
interesting in its own right. Although factorization holds 
beyond doubt in QCD, it is far from obvious that any 
of the regularization schemes used to make a low-energy 
model well defined is compatible with factorization. In 
addition, we want to determine what is the low-energy 
scale, Qo, the model corresponds to. Here we obtain it 
with help of the analysis of PDA and compare it to the 
Qo obtained in deep inelastic scattering (DIS) from the 
corresponding parton distribution function of the pion 
(PDF). 

To a large extent our treatment of PDA parallels the 
calculation of PDF carried out in previous works ]37| , 
|38| , |39|| . There, it has been argued that for inclusive 
processes, such as in deep inelastic scattering, by far 
the most convenient regularization scheme is the Pauli- 
Villars (PV) method. Such a regularization allows the ex- 
traction of the leading-twist contribution to the forward 
virtual Compton amplitude which possesses proper sup- 
port and normalization. The relevance of regularization 
in chiral quark models should not be underestimated; it 
is not evident what is the most convenient way to cut-off 
high energies in such a way that most features of QCD 
are retained. Those include chiral symmetry, gauge in- 
variance, and scaling properties. The main outcome of 
the calculation presented in Ref. ]37j was that, at the 
scale Qo at which the model is defined, the valence PDF 
is a constant equal to one, 

q(x, Qo) = q(l - x, Qo) = V*(x, Q )/2 = 1. (1) 

After QCD evolution at leading order (LO), impressive 
agreement with the analysis of Ref. at the reference 
scale Q = 2GeV has been achieved. At this scale the va- 
lence quarks carry 47% of the total momentum. This im- 
plies a rather low scale Qo, as suggested by the evolution 
ratio a(2GeV)/a(Qo) = 0.15 relevant at leading order. 
For a(2GeV) = 0.32 listed in the PDG @, and for the 
evolution with three flavors, this corresponds to [p"o| 



Qo = 313 MeV, a(Q ) = 2.14 



(2) 



(see Ref. |}7| for details). The low scales are confirmed by 
the next-to-leading (NLO) analysis of Ref. with the 
NLO effects small compared to the LO ones] 51J . Moti- 
vated by this success, in the present paper we investigate 
whether the evolution ratio and the values (||) found in 
deep inelastic scattering are compatible with the values 
extracted from a similar analysis of PDA at LO in the 
same model (NJL) with the same (PV) regularization. 
This is the main objective of this work. 

In the NJL model PDA has already been estimated by 
several authors g ||, |§. The work 01 Refs - H 
uses the Brodsky-Lepage cut-off regularization as sug- 
gested by the light-front quantization formalism. As a 



consequence, the asymptotic form ip(x,Qo) = 6x(l — x) 
is obtained without any additional evolution. On the 
other hand, the same regularization yields PDF of the 
form xV n (x,Qo) ~ 6x 2 (l — x) ^] which is far from 
the asymptotic value xV^{x,oo) = x5(x) = 0. This is a 
rather puzzling result, which may have to do with sub- 
tleties of introducing a regularization in the light-cone 
quantization method (see also Ref. |Q). For that reason 
we prefer to use a manifestly covariant formalism, where 
chiral symmetry can be easily implemented in presence of 
the regularization. In Ref. j2^] PDA has been extracted 
from the transition form factor by examining the asymp- 
totic behavior for large photon virtualitics. This requires 
introducing a regularization for an abnormal parity pro- 
cess which also modifies the chiral anomaly, and hence, 
for typical parameter values f43|] , the 7r° — > 77 decay 
rate is reduced by 40% of the current algebra value. Our 
approach is free of such problems. 



II. THE NAMBU-JONA-LASINIO MODEL 

For the reader's convenience we briefly review the NJL 
model in such a way that our results can be easily stated. 
The SU(2) NJL Lagrangian in the Minkowski space is 
given by @ H 

£ N JL = q(i@ - M )q + ^ ((qq) 2 + {qri^qf) (3) 

where q = (u, d) represents a quark spinor with N c col- 
ors, f are the Pauli isospin matrices, M stands for the 
current quark mass, and G is the coupling constant. In 
the limiting case of the vanishing Mo the action is in- 
variant under the global SU(2)r (gi SU(2)l transforma- 
tions. With help of bosonization, the vacuum-to-vacuum 
transition amplitude in presence of external vector and 
axial-vector currents, (w, a), can be written as the path 
integral 

(0|Texp{i J d^xlqif + ^q])^) 

= J Z?SL>ncxp{i > S'}. 

The following Dirac operators 

?:D = i0 -M - (S + i7 5 T-n)+^ + ^75, 
iD 5 = -i@ - Mo - (S - i-f 5 r- ft) + i> - 4l5, 

are introduced. The fields (S,n) are dynamical, internal 
bosonic scalar-isoscalar and pseudoscalar-isovector fields, 
which after suitable renormalization can be interpreted 
as the physical a and pion fields. The PV-rcgularized 
normal parity (75-even) contribution to the effective ac- 
tion is I35I t5ql 



iN r 



^QtrlogpDg+Af 



it) 



2G I ^ 



n 2 ), 



(4) 
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with tr denoting the trace in the Dirac and isospin space. 
In general, we assume n PV subtractions, with the con- 
ditions J27=o c »A? fe = for k = 0, ...,n , and with cq = 1, 
Aq = 0. At least two subtractions (n — 2), which is the 
case used throughout this paper, are needed to regularize 
the quadratic divergence. The abnormal parity (75-odd) 
contribution to the effective action is 



So 



dd 



-y{trlog(D 2 ) 



trlog(D 2 )} 



(5) 



Notice that no explicit finite cut-off regularization is in- 
troduced in the abnormal parity contribution, as de- 
manded by a proper reproduction of the QCD chiral 
anomaly. This subtle and important point has been dis- 
cussed in detail in Ref. flfjj . 

Any mesonic correlation function can be obtained from 
this gauge-invariantly regularized effective action by a 
suitable functional differentiation with respect to the rel- 
evant external fields. In practice, one usually works in the 
formal limit large N c , in other words, at the one-quark- 
loop level. To fix the parameters in the PV-regularized 
NJL model we proceed as usual (see, e.g., Ref. Q). The 
effective potential leads to dynamical chiral symmetry 
breaking, thereby yielding a dynamical quark mass, M, 
and condensates given by 

(uu) = (dd) = - M 2G M ° = ±N C MI 2 , (6) 

where the quadratically-divergent integral, I2, is defined 

as 



h = i 



d 4 k 



E 



1 



(2tt) 4 ^ 1 {-k 2 + M 2 + Af - ie) 



The calculation of the relevant correlation function yields 
for the pion mass 



2/2 



Mo 



F{ml) M - M 



(8) 



The pion weak-decay constant, and the pion-quark 
coupling constant, g-wqq, are given by 



Trqq, 



1 



n 2 

'J TV- 



AN, 



<i<i 



(9) 
(10) 



respectively. We have introduced the following short- 
hand notation: 



F(p 2 ) = / dxF(p 2 ,x), 
Jo 



(11) 



where, in terms of the PV-regularized one-loop integrals, 



(12) 



1 [-k 2 - x(l - x)p 2 + M 2 + Af - ie} 2 

= (4^)2 E c * lo s l M2 + A * - < l - x )p 2 ] ■ 

The function F in an obvious manner satisfies the sym- 
metry relation F(p 2 ,x) = F(p 2 , 1 — x). In the case of two 
subtractions, and in the limit Ai — ► A2 = A used in this 
paper, we have £. ( cJ(Af) = /(0)-/(A 2 )+A 2 /'(A 2 ). In 
the numerical analysis of this paper we work in the strict 
chiral limit, with Mq = 0. The parameters are fixed as 
usual; we adjust the cut-off, A, in order to reproduce the 
physical pion weak-decay constant, / w = 93.3 MeV. The 
coupling constant, G, is traded for the constituent quark 
mass, M, which remains the only free parameter of the 
model. In our study of the pion light-cone wave function 
we use two sets, which cover the range used in other phe- 
nomenological applications of the model: M = 280 MeV, 
A = 871 MeV (case of Ref. §), and M = 350 MeV, 
A = 770 MeV. These give the quark condensate equal to 
(uu + dd) = -(290MeV) 3 and -(271MeV) 3 , respectively. 
As we shall see, the results are insensitive to the choice 
of parameters. 



III. PION LIGHT-CONE WAVE FUNCTION 
AND PION DISTRIBUTION AMPLITUDE 

The pion light-cone wave function (the axial-vector 
component) is defined as the low-energy matrix ele- 
ment SI 



iy/2 



fit, p i{2x-l)£-p+-Sx-kx 



{TT+(p)\u(C,t±h + l5d(0)\0) 



(13) 



(p n (x) — / d k±^Tr(x,k±) 



where p^ = and p± = 0. The pion distribution 
amplitude is defined as 



(14) 



Formally, in the momentum space, Eq. (|13|) corresponds 
to integration over the quark momenta in the loop inte- 
gral used in the evaluation of f w , but with k + — p + x = 
m^x and k± fixed. Thus, with the PV method and after 
working out the Dirac traces, we have to compute 

iT r 2iN c Mg wgq fdk+dk- 



(2tt)' 



5 (k + — xp + ) 
m v x[\ — x) 



E c ^ x 



(15) 



1 



k 2 + Af 2 +A 2 +iO+ , k 2 +M 2 +A 2 +i0+ 

— - k - 
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where the location of the poles in the k~ variable has 
been explicitly displayed. Evaluating the k~ integral 
gives the pion LC wave function in the NJL model with 
the PV regularization: 



4N c Mg* q 



2^ x 



k\ + A 2 + M 2 - x(l - x)ml ' 



The function is properly normalized, 

/ d 2 k±dx^T T (x,k±) = 1, 

and satisfies the crossing relation 

V n (x,k±) = ^(l-x,kj_). 



(16) 



(17) 



(18) 



For m n 7^ it is non-factorizable in the k± and x vari- 
ables. Integrating with respect to k± yields the pion 
distribution amplitude, 



tp«(x) = AN c MF(ml,x) 



(19) 



The crossing property, (p v (x) = (p„(l — x) follows triv- 
ially, and Eq. (g) gives the correct normalization, namely 
/ dx(p 7V (x) = 1. 

As a consequence of the PV condition with two sub- 
tractions one has, for large 



%ir(x,k±) 



AN C M 2 £V Cl Af 
16tt 3 / 2 k\ 



(20) 



which gives a finite normalization and a finite second 
transverse moment, 



(*5 



d 2 k ± I dx* 7r (a;,fc ± )fci (21) 
o 



In the chiral limit, m n = 0, one can use the Goldberger- 
Treiman relation for the constituent quarks, g^qq]^ — 
M. Then f 2 = 4N C M 2 F(0), which gives the very simple 
formulas 



T . , . 4iV c M 2 v ^ 

**(x,k±) = — — > < 



1 



16^ 3 /2 ^ + A 2 + M 2 



(kl) 



1, 



M (uu) 

J TV 



(22) 

(23) 
(24) 



In the chiral limit ^> 7r (x,k±) becomes trivially factoriz- 
able, since it is independent of x. A remarkable feature is 
that the last two relations, Eq. ( p3| ) and Eq. (|24|), are in- 
dependent of the PV regulators. A similar situation has 
also been encountered when computing PDF in the chiral 
limit ]3j] |; it was a constant equal to one, regardless on 
the details of the PV regulator. We will show below that 
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FIG. 1: The pion light-cone wave function in the chiral limit, 
evaluated in the Nambu-Jona-Lasinio with the Pauli-Villars 
regularization with two subtractions and with the constituent 
quark mass M = 280MeV (solid line) and 350 MeV (dashed 
line), plotted as a function of the transverse momentum k±. 
The wave funcion does not depend on x. The normaliza- 
tion is such that J d 2 k±'^/ n (x, k±) = <Ptt(x) = 1- The second 
transverse moment is (A?5_) = —M(uu)/f% — (625MeV) 2 for 
M = 280MeV and (634MeV) 2 for M = 350MeV . The scale 
relevant for the calculation, as inferred from the QCD evolu- 
tion ||, is Q = 313 MeV. 



by putting together Eq. (|23j) and the results of Ref. |37[] 
an interesting relation follows. 

Higher transverse moments diverge if one restricts the 
number of Pauli-Villars subtractions to two, but Eq. ( |2^ ) 
and Eq. (^4|) remain still valid if more subtractions are 
considered. 

In Fig. (|l]) we show the fc^-dependence of the light- 
cone pion wave function in the chiral limit (finite pion 
mass corrections turn out to be tiny, at the level of a 
few %) for the PV regularization with two subtractions, 
and with M = 380MeV and 350JvleV. For these values 
we get the transverse moment (k 2 ^) — (625MeV) 2 , and 
(634MeV) 2 , respectively. This value is about a factor 
of two larger than the one found in Ref. J25|, namely 
(430MeV) 2 , and a factor of four higher than the findings 
of Ref. (316MeV) 2 , at the scale at which a/ir - 0.1, 
i.e. Q ~ 1 — 2GeV. As we shall see below, a part of 
the discrepancy can be attributed to the QCD radiative 
corrections. 

In non-local versions of the chiral quark model, where 
a momentum-dependent mass function is introduced as 
a physically motivated regulator, the trend to produce 
a constant PDA has also been observed if the constant 
mass limit is considered [^O], ^lj |2f|. In those models 
such a limit effectively corresponds to removing the reg- 
ulator, against the original spirit of the model. Unfortu- 
nately, for the genuine non-local case those calculations 
violate proper normalization of PDA, because the em- 
ployed currents do not comply with the necessary Ward 
identities required by chiral symmetry. The problem has 
been addressed in Ref. p6| , where it has been found 
that about a third of the normalized PDA comes from 
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the non-local currents. For a Gaussian mass function 
there is a clear flattening of <p n (x) in the central region 
of 0.2 < x < 0.8 |||. 

We stress that our result, Eq. (p3|), holds true without 
removing the Pauli-Villars regulator and is in harmony 
with chiral symmetry, since the starting point was the 
normal parity action, which by construction preserves 
chiral symmetry. Obviously, the fact that our final an- 
swer does not depend on the form of the PV regulators 
used makes any subsequent manipulation with the regu- 
lators fully irrelevant. 

Another point is that PDA from Eq. (||) and PDF 
from Eq. (Q) yield the relation (p^{x) — V Tr (x)/2 valid at 
a low scale Qo- It is noteworthy that in the framework of 
QCD sum rules the same identity between PDA and PDF 
has also been obtained |11 at some scale, although there 
the asymptotic form for PDA was assumed without the 
QCD evolution, while PDF was obtained by QCD evolu- 
tion. We will show below that if evolution is undertaken 
for both PDA and PDF at the same low energy scale, an 
overall consistent picture arises. 



O 



9- 



r=0.09 
r=0.21 

6x(1 -x) (Asymptotic Result) 
9^ (x, Q ) = 1 (Initial condition) 




FIG. 2: The pion distribution amplitude in the chiral limit 
evolved to the scale Q 2 = (2.4GeV) 2 . The two values for 
the evolution ratio r = a(Q)/a(Qo) reflect the uncertainties 
in the values of Ref. || based on an analysis of the CLEO 
data. We also show the unvolved PDA, ^^(xjQq) = 1, and 
the asymptotic PDA, ip^ix, oo) = 6:r(l — x). 



IV. QCD EVOLUTION 

The comparison of the leading-twist PDA to high- 
energy experimental data requires, like for PDF, the in- 
clusion of radiative logarithmic corrections through the 
QCD evolution |, |]. For the pion distribution ampli- 
tude this is done in terms of the Gegenbauer polynomials, 
by interpreting our low-energy model result as the initial 
condition. For clarity we work in the chiral limit, hence 

^{x,Q ) = l. (25) 

Then, the LO-evolved distribution amplitude reads [|[ |j| 



n {x,Q) = 6x(l~x)J2C 3 J 2 (2x-l)a n (Q), (26) 



71=0 



where the prime indicates summation over even values of 
n only. The matrix elements, a„(Q), are the Gegenbauer 
moments given by 



a n {Q) 



3/2 - , 

with C n denoting the Gegenbauer polynomials, and 

2 " +1 
" (n+ l)(n±2) ~ 4 ^ k 



2n + 3 



<Q) 

3(n + l)(n + 2) WQo) 



7i°V(2/3o) 



dxCl' 2 {2x- l)<p„(x,Q ), 



(27) 



7< 0) = -2C F 
11 



k=l 



00 



Y C A - -N F , 



(28) 



with Ca — 3, Cf =4/3, and Np being the number of 
active flavors, which we take equal to three [^3|. With 



our constant amplitude ( |25| ) we get immediately 



dxCl /2 {2x- l)(p„(x,Q ) = 1. 



(29) 



Thus, for a given value of Q we may predict PDA. 
We need, however, to know what the initial scale Qo 
is, or, equivalently, to know the evolution ratio r — 
a(Q)/a(Qo). The fitting procedure of Ref. (|] yields 
a 2 (2.4GeV) = 0.12 ± 0.03 (with the assumption a k = 0, 
k > 2). We treat this as experimental input, and then 
with help of Eqs. ( p7|j29[ ) we get for the evolution ratio 



a(Q = 2.4GeV)/a(Q n ) = 0.15 ±0.06. 



(30) 



which at LO implies Qo = 322 ± 45MeV, a value com- 
patible within errors with (Q). 

The fit of Ref. || with non-zero 04 yields a 2 = 
0.19 ± 0.04 ± 0.09 and a 4 = -0.14 ± 0.03 =p 0.09. The 
central value of a 2 would imply, according to our pre- 
scription, the evolution ratio of 0.31, and, correspond- 
ingly, Qo = 0.47to'ig GeV, a much larger central value 
than (H), but with very large errors. For that reason, in 
the numerical studies below we use the value ([50|) for the 
evolution ratio. 

We can now predict the following lowest-order coeffi- 
cients: 



a 4 (2.4GeV) = 0.044 ±0.016 
a 6 (2.4GeV) = 0.023 ± 0.010 
a 8 (2.4GcV) = 0.014 ±0.006 
ai (2.4GeV) = 0.009 ± 0.005 



(31) 
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For the sum of the Gegenbauer coefficients we get the 
estimate 

f>( = 2.4GeV) = f V-feQ^V) ^ 

J 6x(l - X) 

= 0.25 ±0.10 (32) 

where the uncertainties correspond to the uncertainties 
in Eq. ©. 

The leading-twist contribution to the pion transition 
form factor is, at the LO in the QCD evolution [0, equal 
to 



Q 2 F, 



r(Q) 



twist — 2 



\x^ Q \ 
6x(l — x) 



(33) 



The experimental value obtained in CLEO [|| for the 
full form factor is Q 2 F r ^ 1 (Q)/(2f 7r ) = 0.83 ± 0.12 at 
Q 2 = (2.4GeV) 2 . Our value for the integral, 1.25 ± 0.10, 
overestimates the experimental result, although at the 
2cr-confidence level both numbers are compatible. Tak- 
ing into account the fact that we have not included nei- 
ther NLO effects nor an estimate of higher-twist contri- 
butions, the result is quite encouraging. 

In Fig. |2| we show our PDA evolved to Q — 2.4GeV, 
for two values of the evolution ratio, which reflect the 
uncertainties from Eq. (30). We also show the initial 
and the asymptotic PDA's. It is interesting to note that 
after evolution our results closely resemble those found in 
transverse lattice approaches [jl6], [l7], In particular, 
we get for the second ^-moment (£ = 2x — 1), 

(C 2 ) = / dx(p„(x,Q = 2.4GeV)(2x- l) 2 
Jo 

= 0.040 ±0.005, (34) 

to be compared with (£ 2 ) = 0.06 ± 0.02 obtained in the 
standard lattice QCD for Q = l/a = 2.6 ± O.lGeV 
From the PDF calculation at LO of Ref. j37j we estimate 
that if the momentum fraction carried by the valence 
quarks at Q = 2GeV is 0.47 ± 0.02%, then Q is such 
that a(Qo) = 2.14, and the evolution ratio at Q — 2GeV 
is r = 0.15. Then, for Q = 2.4GeV we get r = 0.14 from 
the analysis of PDF, a value compatible, within uncer- 
tainties, with the present calculation, Eq. (|30|). This is 
a crucial finding, showing the consistency of the results 
obtained in our approach. 

One might worry that the starting condition ( |25| ) 
does not satisfy the end-point vanishing behavior and 
therefore cannot be expanded in terms of the Gegen- 
bauer polynomials. This is true, provided one insists 
on uniform pointwise convergence. However, the Gegen- 
bauer polynomials form a complete set in the space of 
square-summable functions, hence convergence may be 
understood in a weak sense [Q. The slow convergence 
is reflected by the fact that in Fig. || at least 30-100 
Gegenbauer polynomials are needed for evolution ratios 
r = 0.9 — 0.21 respectively. The convergence at the mid- 
point, x — 1/2, is improved, since the series for tp(x,Q) 



is sign-alternating. At the end-points, x = 0,1, the se- 
ries diverges, since Cl /2 (±l) = ±(2fc + l)(2fc + 2), which 
means that the convergence in Eq. ( |26| ) is not uniform. 
In order to analyze the behavior close to the end-points 
in a greater detail we consider the large-n contribution 
to Eq. ©. We have 



<*(Q) 
a(Qo) 



7i 0) /(2/3o) 



3d 



In - 



(<>()) . 



(35) 



hence, for Q — ► Q , Q > Q , and with x — > (recall that 
the function is symmetric under x — > 1 — x), we obtain 



ip^ix -> 0,Q) 



V 2/3o a(Q ) 



(36) 



71=1 
1 _ 



i z is the Riemann ( function, and 
oo. Thus the slope of the evolved 



where C(z) = 

C(i) = EZ ~ 

PDA at the end-points becomes steeper and steeper as 
Q -> Oo- 

The QCD evolution also influences the value of the 
transverse moment. According to the work of Ref. [ E5[ , 
{k]_ ) can be expressed as (A: 2 ) = 5mQ/36, where m 2 , = 
{qa^ v F^ v q) / (qq) is the ratio between the quark-gluon 
and quark condensates. The quantity m 2 , is scale de- 
pendent and has been estimated to be mg(lGeV) = 
0.8 ± 0.2GeV 2 [|6). Using the corresponding anoma- 
lous dimensions, 4 for (qq) and —2/3 for (qa^ F^q) p7| , 
yields 



(kl) Q = (*(Q)\ (4+2/3W ° 
(ki) Qo \a(Q )J 

a(Q) \ ^(33-2^) 

a(Q ) 



(37) 



For Nf — 3 this scale dependence can be seen in 
Fig. (|). For the values Q = 1 - 2 GeV one gets a 
reduction factor of 0.37 — 0.45 for the ratio (|37j), and 



(k 2 ± ) Q = (430MeV) 2 - (380MeV) 2 for the second trans- 
verse moment, somewhat higher than the QCD sum rules 
estimate based on Ref. @, (316MeV) 2 , or on Ref. ^6), 
(333±40MeV) 2 . 



V. THE RELATION TO DEEP INELASTIC 
SCATTERING 

As we have already stated in Eq. (|l|) , the valence PDF 
for the pion in the chiral limit has also been found to 
be a constant equal to one p7j . At LO the non-singlet 
evolution of the PDF moments is quite similar to that of 
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FIG. 3: Dependence of the second transverse moment of 
the pion light-cone wave function, (fcj_)<3 / {k\)Q Q (solid line), 
the second Gegenbauer moment a 2 (Q) /a2(Qo) of the pion 
distribution amplitude (dashed line), and the evolution ra- 
tio r = a(Q)/a(Qo) (dotted line), plotted as functions of 
the scale Q. The leading-order QCD evolution is applied. 
All quantities are relative to their values at the low energy 
scale, Qo = 313 MeV, at which the momentum fraction car- 
ried the quarks equals unity |}7|, according to the prescrip- 
tion that in a quark model Qo is denned by the condition 
(iKr(x, Qo)) = 1. In our model a(Q ) = 2.14, a 2 (Qo) = 7/18, 
and (fei)Q = (625MeV) 2 for M = 280MeV and (634MeV) 2 
for M = 350MeV and in the chiral limit. 



the Gegenbauer moments of PDA, Eq. (|2q ) , namely 

a(Q) 

2 ( a(Q) 



(38) 



dx x n V 1T (x, Q a ) — 

7i°V(2/3o) 



Thus, for n = 2, one obtains 
a 2 (Q) (x 2 V w (x,Q)) 



a.2(Qo) {x 2 V^{x,Qq)) \a{Qo) 



7fV(2/3o) 



(39) 



For Np = 3 this scale dependence for the ratios can be 
looked up in Fig. (§). Using (x 2 V 7T (x,Q )) = 2/3 and 
a 2 (<9o) = 7/18 yields 



02 (Q) 



(x 2 14(x,Q)) 



7_ 

12' 



(40) 



hence a 2 (2GeV) = 0.12±0.01 for (x 2 V n ) = 0.20±0.01 M 
and a 2 (2GeV) = 0.10±0.01 for (x 2 ^) = 0.17±0.01 @|. 

One can combine Eqs. (0,^,p6|,^9|) to obtain the fol- 
lowing very interesting LO relation that holds in the con- 
sidered model: 



<Ptv(x,Q) 
6x(l — x) 



dyK(x,y)Vir{y,Q), 



(41) 



where the kernel K is independent of Q 2 , and is given by 



K(x,y) 



E 3(ra 

n=2 y 



(2n + 3) 
~2) 



Cl' 2 {2x - l)y r - 



(42) 



In general, the relation (|4l]) holds in any model where 
PDA and PDF are simultaneously equal to unity at some 
scale Qo, and are subsequently evolved at LO. Physi- 
cally, Eq. ( f4l|) simply tells us that the departure of PDA 
at a given Q 2 from the asymptotic form is proportional 
to a weighted integral of PDF at the same Q. Clearly, 
tp w {x,Q) — > 6a;(l — x) if V w (x, Q) — > 2S(x) or, equiva- 
lently, xV 7T (x, Q) — > 0, since K(x, 0) = 0. Roughly speak- 
ing, in the present model the pion distribution function 
is as close to the asymptotic value as the non-singlet par- 
ton distribution. A remarkable feature of relation (JiTJ) is 
that it binds matrix elements related to exclusive (PDA) 
and to inclusive (PDF) processes. 

In order to evaluate the kernel we use the symmetrized 
generating function of the Gegenbauer polynomials, 



G(x,y) 



E 

n=2 



Cll 2 {2x~lW> = \{R-^ + R^} 



- 1 



R± = lT2(2x-l)y + t/ 2 : 

whence one can obtain 

2 1 

K(x,y) = -G(x,y) - 

The integrals can be worked out to yield the final result 

1 



(43) 



dy'y'G(x,y'). (44) 



K{x,y) 



24i?l / V(z 



[8 (x-l)xy 2 + R+((2x-l)y-l) 

+2 v /^7 (a - 1) xy 2 (1 + (2 - 4a?) y + y 2 ) 

+Rl(l-8(x-l)xy 2 )]- (y~-y) (45) 

To test the success of Eq. ( fi"l| ) we need some input for 
V v (x,Q). However, taking into account the fact that the 
agreement of the evolved valence PDF, V !V (x, Q) with the 
parameterization of Ref. |4(| at Q 2 — 4GeV 2 is almost 
perfect |37], ^9|, and that the results are almost insen- 
sitive to the evolution ratio, a(Q)/a(Qg), Fig. || can be 
regarded as a direct prediction of Eq. (|l]) taking Ref. p(J 
as input for V^(x, Q). A further consequence of Eq. pi] ) 
may be obtained by integrating with respect to x and 
performing the sum over n. Through the use of Eq. ([32] ) 
we get 



oo ' -i 

V] a n (Q) = / dyK(y)V n (y,Q) 

n=2 J ° 



where 



«(y) 



„1 oo ' 

/ dxK(x, y) = V 

J0 „-o 



(2n + 3) 

^ 2 WT2) 



(46) 



y n (47) 



3y 2 + l log(l - y) + log(l + y) 



6(l-y 2 ) 



6y 2 
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Notice that, for Q — > oo we get V K (x,Q) — ► 26 (x) and 
since n(y) = 7y 2 /12 + 0(y 4 ) one gets Yln=2 a n(Q) ~* 
0, as expected. Finally, using the parameterization of 
Ref. (40) wc get @ 

OO ' 

a n (2GeV) = 0.25 ± 0.03, (48) 

a value perfectly compatible with Eq. (|32| ) although with 
smaller uncertainties [^6| . Again, this verifies the consis- 
tency of our approach. 

VI. CONCLUSIONS 

We summarize our points. We have computed the 
light-cone pion wave function and the pion distribution 
amplitude in the Nambu-Jona-Lasinio model. To this 
end, and to comply with previous results regarding the 
parton distribution functions, we have used the Pauli- 
Villars regularization method in such a way that chi- 
ral symmetry, gauge invariance, and relativistic invari- 
ance are preserved. As a result, we find that in the 
chiral limit the pion distribution amplitude, computed 
as a low energy matrix element of an appropriate op- 
erator, is a constant equal to one, </? w (x) = 1, and the 
second transverse moment of the pion light-cone wave 
function is {k\) = —M(uu)/f%, with M denoting the 
constituent quark mass. Both results are independent of 
the particular form of the Pauli-Villars regulators used. 
After the QCD evolution of the pion distribution am- 



plitude to the experimentally accessible region we find 
a result still rather far away from the asymptotic form, 
ip n (x) = 6x(l — x), but in a good agreement with the 
analysis of the experimental data from the CLEO collab- 
oration. We can determine the working momentum scale 
for the model to be Qo = 313 MeV, a rather low value. 
Moreover, the scale Qo obtained in this work is com- 
patible, within experimental uncertainties, to the value 
obtained from the previous analysis of the parton dis- 
tribution functions, carried out within exactly the same 
model. At the scale Qo the quarks carry all the mo- 
mentum of the pion. Our value obtained for the second 
transverse moment of the pion light-cone wave function, 
(kj_), becomes, after the QCD evolution, not far from the 
estimates based on the QCD sum rules. Finally, we have 
also derived a model relation which binds the departure 
of the pion distribution amplitude from its asymptotic 
value to an integral involving the pion quark distribu- 
tion function. The relation, specific to the feature of 
our model that at the scale Qo both the PDA and PDF 
are constant and equal to unity, has been successfully 
checked against the available data. 
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